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Abstract 

We present quantum stochastic calculus in terms of diagrams taking 
weights in the algebra of observables of some quantum system. In par- 
ticular, we note the absence of non-time-consecutive Goldstone diagrams. 
We review recent results in Markovian limits in these terms. 
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1 Introduction 

Quantum stochastic calculus (QSC) involves an analysis of the fundamental 
quantum processes of creation/annihilation/conservation 1. and intuitively this 
is somehow related to emission/absorption/scattering of physical quanta as de- 
scribed by quantum field theory (QFT). Quantum stochastic theory has the 
advantages, as well as the limitations, of having a mathematically rigorous set- 
ting. It also has the theory of classical probability to fall back on for much of its 
inspiration. So much so that the relationship with QFT, which was originally 
a major motivating factor, is now frequently overlooked. Effectively, the fun- 
damental quantum processes should be idealizations of quantum fields for some 
suitable "Markovian" regime. They were introduced to describe open systems 
dynamics: here the quantum noise couples to some quantum system and so, in 
some sense, we are dealing not just with traditional quantum fields, but with 
quantum fields taking values in the algebra of observables of some quantum 
system. 

When presenting his famous list of problems, Hilbert is supposed to have 
quoted an unnamed colleague as saying that "a mathematical theory should 
not be considered complete until one can walk out into the street and explain 
it to the first person you meet" . Let us suppose we did this and, as luck would 
have it, the first person we meet is a physicist. Would we succeed in explaining 
quantum stochastic calculus? My contention is that we should, though we might 
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have to make do with some formal mathematics (presumably Hilbert wouldn't 
have objected?). There are many fundamental ideas, familiar to physicists, 
hidden (sometimes too well-hidden!) in the mathematical formalism of quan- 
tum stochastic calculus. In this article, I've tried to present some basic results 
from quantum stochastic calculus in the language of QFT and, in particular, 
in terms Feynman-type diagrams: the hope being that mathematicians and 
physicists will learn something from the cross-over. Diagrammatic techniques 
are central to QFT yet also of independent interest from it; they still re- 

main an essential tool for investigating mathematical aspects of quantum theory 
and continue to yield some of the most illuminating insights (see, for instance, 
the papers ^ El)- I also want to present an account of a recent paper [7] 
which deals with the QSC approximation and which was originally formulated 
in diagrammatic language. 



1.1 Expansions of Evolutions 

A free-particle (Bosonic) quantum field $ t living on a Hilbert (Fock) space S) can 
be decomposed into positive and negative frequency terms as $ t = $j ^: 
here we suppress all dependence other than time and understand the time label 
to refer to Heisenberg picture of 
annihilation field and ^[~^ to be 
vector, then we have the identity 



to refer to Heisenberg picture of the free dynamics. We take 3>j to be an 
annihilation field and $^~' ) to be a creation field. If Cl is the Fock vacuum 



H = 0, (1) 
along with the canonical commutation relations 

G(t,s). (2) 



Here G (t, s) = (f2|$t$ s CI). Related to this is the propagator K defined as 
K (t, s) = Ul\f $ t $ s ft) = G (t, s)9(t~s) + G (s, t)9{s-t). 



As usual f is Dyson's chronological operation placing Heisenberg picture oper- 
ators in increasing time-order from right to left. 

Let {T t : t > 0} be a family of self-adjoint operators on Sj with T f being 
some function of <5j We are then interested in the evolution operator 

U t =f jexp-i / T s ds\ , (3) 
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by which we mean the solution to the Schrodinger equation idtVt — TT t £/t, 
Uq — 1. Suppose that we have a polynomial interaction 



2 



The standard device of quantum field theory is to expand U t in terms of dia- 
grams, see, e.g. |2] or A Wick diagram D is constructed as follows: choose 
n = n(D) labelled points (vertices), each vertex will have some labelled legs 
attached (we let m v = m v (D) denote the number of vertices having v legs so 
that n — J2u m v)i we j° m several pairs of legs to form (undirected) edges, the 
result is a graph having several external lines and we now ignore the labelling. 
A Wick diagram is then the class of all topological equivalent graphs. We let 
c = c (£>) denote the number of ways we could have originally connected the 
various legs to get the same graph. Now denote by T>w the set of all Wick 
diagrams and define, for each D 6 T>w, the operator 

6 w = h>" £ n r * Lr 9 K n (* (+) + $( i (4) 

where N is normal ordering (placing all creation fields to the left of all 

annihilation fields $( - )) and under the integral we have a factor K (ti,tj) for 
each edge (i,j) occurring and a factor $ tfc for each external line at vertex (k). 
It is then a basic result of QFT that Ut admits the expansion 

u t = E ^w- 

z?ex> w 

Next let denote the subset of connected Wick diagrams then we may list 
the elements as Pi , P%, ■ ■ ■ and each D 6 T>w can be decomposed as D = 

P (t) ni P (t) n2 

P™ 1 x P™ 2 x • ■ ■ . Now one readily checks that D (t) = ^ 

ni! n 2 ! 

and so 

ni,n 2 ,"'i=0 J j=0n=0 J P£?w 

What we have managed to do is to express the evolution operator Ut as a normal 
ordered exponential of a sum over connected Wick diagrams. The connected 
Wick diagrams play the role of the 'primes' amongst the set of all Wick diagrams 
- indeed the trick of replacing a sum of products by a product of sums is just the 
one that goes on when we develop a prime number expansion of the Dirichlct 
series of a multiplicative function, the Riemann zeta function being perhaps the 
best known example; it is also the trick used to compute the grand canonical 
partition function for the free Bose gas. The result should be understood as an 
operator theoretic version of the usual cumulant moment expansion. 

Now let t)s be a fixed Hilbert space. We move the action up to the Hilbert 
space f) ® S) and set 

r f = E^KTH +) f ^ 

a, 13 
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where we take E^p = Ep a . We now introduce a class of diagrams known as 
Goldstone diagrams - they differ from the previous ones in that the vertices are 
placed in time order ,2 . Consider times t n > ■ ■ ■ > £2 > tl in the interval [0, t] 
and draw these as vertices as shown below: 



Suppose that at vertex j we have (3j legs coming in from the right, repre- 
senting annihilators, and atj legs going out to the left, representing creators. For 



example, the vertex for E23 



(-) 



and [3 ■ = 3 annihilators is sketched as 



(+) 



where we have a q ■ = 2 creators 



aj creators 




/? annihilators 



We construct a Goldstone diagram D as follows. We take an arbitrary 
number n — n (D) vertices and draw in an ordered line as above. We then 
draw creation / annihilation lines at each vertex corresponding to one of the 
terms appearing in T. We then connect selected creation legs to (necessarily 
later time) annihilation legs: the remaining uncontracted legs are then directed 
external lines. We consider the family T>q of all (topologically distinct) diagrams 
obtained in this way. To each Goldstone diagram D we associate the operator 



d (t) = Hr s anPn . . . E ai/S± ® / n *(-) n g n 

JA n (t) D D D 



(7) 



where A n (t) is the simplicial region {(*«,••• , ti) : t > t n > ■ ■ ■ > t 2 > t\ > 0} 
and we have a factor G (ti — tj) for each edge note ti > tj, a factor <$>^' 

for each incoming external line to a vertex k and a factor for each outgoing 
external line. We then find the expansion 



U t 



DeV G 



(8) 



In QFT one is used to switching between an expansion in terms of Gold- 
stone diagrams (JBJ and one in terms of Wick diagrams 10). In the present case 
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however we have an obstruction: the E a ^s do not necessarily commute! This 
complication means that the Goldstone diagrams are more fundamental in the 
present case. The problem, of course, is that the Dyson operator T is reordering 
the Hcisenberg fields only, while the E a ps remain in their original order. 

1.2 Zero Dimensional QFT 

Let a and a 1 be annihilation operators for a single mode harmonic oscillator, 
We have the commutation relations [a, aH = 1 and aft = 0. Let us consider the 
observable q = Z(V + z*a where z is a complex number. The Baker-Campbell- 
Hausdorff theorem says that exp {itq} = exp {itza*} exp {-±\ z \ 2 t 2 } exp {itz*a} 
which here has the same content as the expansion JSJl. We can use a diagram- 
matic presentation based on two types of vertex: the creation type 

which has weight z, and the annihilation type ,. 

r. 

which has weight z* . If 

we take a vacuum expectation of exp {itq} then we need only consider connected 
diagrams having no external lines - and there is only the one! The cumulant 
expansion is then 

(n|exp{itg}fi) = exp{-^ {Z^\ } = exp {~\t 2 \z\ 2 }, 

and we see that q is Gaussian in the vacuum state. As is well known, the 
odd moments vanish while the even moments are (p<\q 2n ty — \z\ 2n ^r^ where 
the combinatorial factor counts the number of ways to partition the 2n (time- 
ordered) vertices into n (contraction) pairs. For instance, the fourth moment 
involves three disconnected diagrams (our convention is that we only consider 
the contractions as edges; the thin horizontal base line does not affect connec- 
tivity!) 



(<? 4 ) = ....cn..cn.... + J..CXI....+ ....an.... = 3M 4 . 

We may also consider the variable N = (a + z)^ (a + z) = a^a + za) + z*a + 

\z\ 2 . We now introduce two extra vertices: a scattering vertex .T^y^T. with 

weight unity and a constant vertex with weight \z\ 2 . To determine 

the vacuum expectation of exp {UN}, we once again sum over all connected 
diagrams with no external lines. This gives 



<n| exp {/7.V[ Q) = exp{M + ...,/ZY... +^ 
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The n th term in the exponential will look like ....{Z^C^C^C ' ' '~YZY-. and 
each such term has weight \z\ 2 - since the scattering vertices all have weight 
unity. All cumulants are equal and we therefore have 

(n\exp{itN}Sl) =cx P |^^-|z| 2 | =exp{|z| 2 (e (t -l)} 

and we recognize N as having a Poisson distribution of intensity \z\ 2 . The 
moments of the variable N are given as a polynomial of degree n in \z\ 2 , vis. 

n 

(n\N n n) = S(n,m)\z\ 2m 

m— 1 

and, as is well-known in combinatorial analysis UJ, the coefficients S (n,m) 
are the Stirling numbers of the second kind: they count the number of ways 
to partition n items into m non-empty subsets. To see why they arise here, 
consider the following diagram contributing to (N 7 }: 



. I Y* Y* i"TT i . 

tj t@ k ti ti ti 

This diagram partitions the 7 vertices into 3 subsets, namely {tr, te, t$, £3}, 
{ti,ti} and {£2}, with each subset forming a connected sub-diagram. This 
contributes (\z\ 2 ) to (iV 7 ). Consulting a textbook on combinatorics to get 

the Stirling numbers, we find (N 7 ) = Y7 m =i S ( 7 > m ) (NT" = \ z \ 2 + 63|z| 4 + 
301|z| 6 + 350|z| 8 + 140|z| 10 + 21\z\ 12 + \z\ u . Alternatively, we could draw all 
Br = 877 diagrams out! The numbers B n = X)m=i & ( n i m ) counting the total 
number of ways to partition the n vertices into non-empty subsets (of connected 
Goldstone diagrams) are known as the Bell numbers. 



2 Quantum Stochastic Calculus 

Remarkably, the equivalence between expansions (jSJl and JSJl is restored in the 
non-commutative case in one very important situation. This is when we consider 
the cases a and j3 taking only the values 0, 1 in JHJ and when the two-point 
function G is replaced by a delta-function. Effectively the field is some form 
of quantum white noise in time. As a, (3 is restricted to either or 1, we 
shall have only four types of vertex: a constant vertex Eqq (g> 1, an emission 
vertex E\$ (g> an absorption vertex Eq\ <8> <I>( + ) and a scattering vertex 

En (g) $(")$(+). 

The reason for the algebraic equivalence, despite the fact that the E a p need 
not commute, is that many of the Goldstone diagrams vanish identically. This 
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is because to the singular nature of the two-point function with respect to the 
simplicial integration in Q . We note that absence of certain diagrams describing 
moments of quantum noises has occurred elsewhere, in particular, there is an 
elegant description of the various forms of independent quantum processes in 
these terms [TT] . 

Let us introduce some specific notations |10|. We make the replacements 
at, ^ a\, G (t, s)^g(t-s), 

where g (t — s) = Kd (t — s) + n*d (t — s). Here k is a complex damping constant 
with 7 = 2Re{fi} > 0. The objects D± (t — s) are one-sided delta functions 
defined by / / (s) £)± {t-s)ds = f (i ± ) 

f(s)d ± (s-t)ds= / f(t + u)* ± (u)du = f(t ± ). (9) 

Let us briefly indicate how to convert Ut = Texp (—i J Q * T s ds| to normal order 

[TO] where T t = £7 ai g ® ^a|^j (a*)' 3 (we use a convention from now on that 

repeated Greek indices are summed over values and 1). When evaluating 
Goldstone diagrams, we find that if the contractions are not time-consecutive, 
that is, if we encounter g (tj — tj) with i > j + 1, then we force the multiple 
equalities = ti-\ = ■ ■ ■ = tj+i = tj due to the time ordering, and so the contri- 
bution vanishes. Only Goldstone diagrams with time-consecutive contractions 
are non-zero. 

Starting from the integro-differential equation Ut = 1 — i L T s U s ds, we have 

[a t ,U t ] = -i [a t ,T s ]U s ds 
Jo 

= -i f 9 (t - s) Em (a t f U s 
Jo 

= -inE x $ (a t f U t 
or a t U t = (1 + ikEu)^ 1 [U t a t — inEioUt] and so 

d t U t = -iE aP ® (a])" (a t f U t = (atf L a0 U t (a t f (10) 

where 

= —iE a p - K,E a i— — — —Eip. (11) 

We may interpret the conversion of the Schrodinger equation to normal ordered 
form as a change from a Stratonovich to an Ito description. This agrees with 
the interpretation given originally by von Waldenfels for emission-absorption 
interactions |12) . 
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Having normal-ordered the Schrodinger equation, we now iterate to get 



U t = 1+ / (at) a L af3 U s (a s fds 
Jo 

= E j t K) a "---(<) ai (W„---i^ 1 )(a* 1 ) /3l "--(citJ 



(12) 



This is reasonably familiar to quantum field theorists and such expressions can 
be found for instance in Berezin's book If / is a suitable test function, wc 
may consider its coherent (i.e. exponential) vector |e (/)) and take at \e (/)) = 
f(t)\e(f)) and (e(f)\a\ = (e (f)\ f (t)* . As (fT2l is normal ordered, we have 
no difficulty in assigning a meaning to (e (f) \ Ute(g)). At this stage we could 
just as well take (|12fl as the definition of the process, this is the starting point 
of the Maassen kernel calculus [TJ. As such the time-consecutive contraction 
property is built into QSC, though in a way that is not readily apparent. 

For the benefit of quantum probabilists, who may well be a little lost at this 
stage, we convert (ITTTf) into more familiar language pQ. Let A" 13 = f* (ajj)" (as)' 3 ds 
and we interpret these as the four fundamental quantum processes: A°° is time, 
Aj° is creation, A^ 1 is annihilation and A^ 1 is conservation. Loosely speak- 
ing, we say {X t t > 0} is adapted if [a;, X t ] — whenever s > t. Setting 

X( = Jq ( a t) a x al ( s ) ( a s)^ ds where the x^p (■) are adapted, we see that 
putting to normal order yields 



* t (1) 4 2) = (at 



X^x^l (s) + x« (s) XP + x$ (s) (s)] (a s f ds. 

The basic idea goes back to Hudson and Streater JS]. In QSC, we usually write 
dX t — x a j3 (t) dA"^ and the above result is presented as the quantum Ito formula 
d (X^XW) = XWd (X&) + d (If 1 )) XW + d (XW) d (XW) along with the 
quantum Ito table dA^dA^ = rfA"' 3 . The equation ifTUI) is then interpreted 
as the Ito quantum stochastic differential equation dllt — L a pUtdA"^ with 
Uq = 1. The coefficients satisfy the identities L a/ 3 + L^ a + jL\ a Lip = which 
are necessary and sufficient for Ut to be an adapted, unitary quantum stochastic 
process. The formula for the product of several quantum integrals comes down 
to a normal ordering problem which can ultimately be presented as a sum over 
diagrams, or equivalently, a sum over partitions of the time indices: for the 
classical case, see [TB] . 



3 Markov Limits 

Finally, we wish to comment on how regular quantum fields can approximate 
the singular fields considered above. Let A ^ be a parameter and consider 
fields (A) with a regular two-point function G\ (•) which becomes a delta- 
function in the limit A — > 0. In particular, we may take G (•) to be a integrable 
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function with 7 = G and k = J °° G and assume that G{— t) = G (t)" 



Then set G\ (t, s) = \~ 2 G j . We would the n argue that in the limit G\ 

converges to the singular function q consider above. We consider the regular 
unitary evolution operators 

U t (A) = fexp j-i jf E aP ® (A))" ($W (A)) 

and we claim that for bounded -E a| a, with ||k£^h|| < 1, Ut (A) converges to the 
singular process Ut considered in the last section. 

The first remark that we make is that the A — > limit leads to the vanishing 
of each non-time-consecutive Goldstone diagram contributing to Ut (A). More- 
over, when the surviving terms are computed and re-summed, we formally get 
the correct Ito expansion l|12l) . Note that L a p = —iE a p—i Y2rLi ^ai {—iK>E\xf 1 Eip 
giving the contribution to a time-consecutive block with a outgoing, (3 incom- 
ing lines and a sum over r successive scattering in between. We see that the 
condition ||k2?h|| < 1 is necessary to sum the geometric series. 

The re-summation is rather tedious, though it helps that we know what 
answer to expect! We also have the issue of convergence, however, we settle 
this below. We remark that it is sufficient to consider only the vacuum conver- 
gence as the more general situation can be inferred from this when we look at 
convergence in arbitrary but appropriately scaled coherent vector states. 

3.1 Pule Inequalities (Gaussian) 

Let us start with the case where we have emission and absorption only in the 
interaction. The vacuum Goldstone diagrams, as we have seen in section 2, 
consist of ri2, say, pair contractions only. A typical diagram, one of 
having 2«2 vertices, is sketched below for 712 = 6: 



. *7 : ; *777» i . . * . 

There exists a permutation a of the n = 2n 2 time indices which re-orders to 
the diagram Dq (n) shown below 



ta(n) t<j{2)ta(l) 



9 



The permutation is moreover unique if it has the induced ordering of the 
emission times. Not all permutations arise this way, the ones that do are termed 
admissible. We now consider an estimate of the n— th term in the Dyson series: 



ni^i = e / n \^ 



O fj 



DeT> G JA i( t ) D Admissible permutations A "(*) _D (n) 



/ n \Gx{t2k-t 2k -i) 

J Bit) k=l 



where R (t) is the union of simplices {(t n , ■ ■ ■ ,ti) : t > £ CT -i(„) > • ■ • > t CT -i(i) > 0} 



over all admissible permutations a. R(t) will be a subset of [0,i] 2 ™ 2 and if we 
introduce variables t 2 k and s 2 k = t 2 k — t 2 k-i for k = 1, • • • , n 2 it is easily seen 

max (i, l)" 2 

that the above is majorized by |k|" 2 x . This the Pule inequality 

n 2 \ 

|17| and clearly gives the uniform absolute estimate required to sum the series. 
3.2 Pule Inequalities (Poissonian) 

We now consider scattering, and constant, terms in the interaction [7|. As we 
have seen there will be B n (the n-th Bell number) Goldstone diagrams con- 
tributing to the n— th term in the vacuum Dyson series expansion. The Bell 
numbers grow rapidly and have a complicated asymptotic behaviour. The pro- 
liferation of diagrams is due to mainly to the multiple scattering that now may 
take place. 

Let us consider a typical Goldstone diagram. We shall assume that within 
the diagram there are ri\ singleton vertices [• ■ • • • • ], n 2 contraction pairs 

[• • • ■ ■ ■], 71.3 contraction triples [• • • ■ ■ ■], etc. That is 

the Goldstone diagram has a total of n — Y]. jrij vertices which are partitioned 
into m = Y]j rij connected subdiagrams. For instance, we might have an initial 
segment of a diagram looking like the following: 



; ; ; . r~r\ y \ v 7\ 



There will exist a permutation a of the n vertices which will reorder the 
vertices so that we have the singletons first, then the pair contractions, then the 
triples, etc., so that we obtain a picture of the following type 

ris triples n 2 pairs n\ singletons 
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The permutation is again unique if we retain the induced ordering of the 
first emission times for each connected block. We now wish to find an uniform 
estimate for the n-th term in the Dyson series, we have 

J2 I TT |Ca| x "weights" (13) 

,• JA„(t) 



Goldstone diagrams 



where the weights are the operator norms of various products of the type 
E ari p n ■ ■ ■ E ai p i . In general, the weight is bounded by 

I l-^n ||" 1+2 " 2+3n3+ '" x (7 n i+™2+™3+--- 

where C = max a/ g ||-E a( 3||. This is because each connected diagram of j vertices 
will typically have one emission and one absorption, but j—2 scattering vertices. 
The Pule argument of rearranging the sum over diagrams into a single integral 
over a region R (t) of [0, t] n again applies and by similar reason we arrive at the 
upper bound for (|13fl this time of the type 



iim+2n 2 +3n 3 +- x Qm+ni+n 3 +- x max L ) 



Here the sum is restricted so that ■ jrij — n. An uniform estimate for the 
entire series is then given by removing this restriction: 

exp{j: (Aj+B)nA 

3(A,B) = V , , , >- 

^ ni!n 2 !n 3 !--- 

where e A = and e B = Cmax(t, 1). Again we use the trick to convert a 

sum of products into a product of sums 

s(A b) = V TT exp ^ J + B ^ n ^ =]JJ2 GXP ^ + B ^ ^ 



ncx P { e (w} 



ex P <; J2 eA3eB 



exp • 



j 

e A+B 
l-e A 



where we need e A < 1 to sum the geometric series - this however, is precisely 
our condition that IIk-ExiII < 1. 
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4 Conclusions 



We have established a Markov limit in the sense of JS| which we may write as 




On the left hand side we have an unitary which can be expanded as a nor- 
mal ordered expression of the quantum fields in terms of Goldstone diagrams. 
The right hand side can be developed as an expansion over time-consecutive 
can be understood as Hudson-Parthasarathy unitary quantum stochastic pro- 
cess. We have shown the non-time-consecutive terms on the left hand side make 
a negligible contribution in the Markovian limit. Interpreting Weyl order as 
Stratonovich form and Wick order as ltd form, the above result can be consid- 
ered as a non-commutative version of the Wong-Zakai limit theorem for classical 
processes. 

The same holds for Fermi fields, however, the proof is complicated because 
we have to take the limit in matrix elements of appropriately scaled number 
states 221- The same basic estimates suffice once more and in the limit we end 
up with the same process except with the A Q ^ now being Fermionic noises. As 
one might suspect, we have to bother ourselves collecting factors of —1, and 
one would expect to obtain the same result if we dealt with q-commutation 
relations PH]. 

We remark that the time-ordered exponentials developed in |21| differ from 
the notions presented here, as we are time-ordering quantum white noises and 
not Ito differentials, though they do arise in models for Markov limits of discrete 
time systems |22|. 

Finally, we mention that we also have the convergence of the Heisenberg 
dynamics U t (A) f (X ® 1) U t (A) to u} (X ® 1) U t 0. This requires a slightly 
deeper analysis, however, the basic estimates above are again at the heart of 
things. We invite the reader to try and imagine the Goldstone diagram expan- 
sion of Ut (A) (X ® 1) U t (A) to get an idea of what is involved. 
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